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ABSTRACT 

We present the integrated Ha emission line profile for 157 Hn regions in the cen- 
tral 3.4' x 3.4' of the galaxy M 83 (NGC5236). Using the Fabry-Pcrot interferometer 
GHaFAS, on the 4.2 m William Herschel on La Palma, we show the importance of 
a good characterization of the instrumental response function for the study of line 
profile shapes. 

The luminosity-velocity dispersion relation is also studied, and in the log(L) — 

_i_ . log(a) plane we do not find a linear relation, but an upper envelope with equation 

• ■ log(LHa) = 0.9 •log(er) +38.1. For the adopted distance of 4.5 Mpc, the upper envelope 

*Q | appears at the luminosity L = lfj 38 - 5 ergs" 1 , in full agreement with previous studies 

of other galaxies, reinforcing the idea of using Hn regions as standard candles. 
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1 INTRODUCTION great efforts h ave been made to e xplain the line profiles 

of Hn regions. iMelnickl (1 19771 ) and iTerlevich and Melnickl 
H n regions are among the brightest non-ephemeral objects l lr . 01 k , 110 . , , . , , „ , ,,. 

° , , (I1981I) analyzed 12 regions and claimed to find a relation 

in galaxies, with prominent emission lines characterizing , ; . , . ., . , , . . , , .... 

between sizes, luminosities L, and integrated emission line 
their spectra. By studying the integrated line prohles of , ., ,. . fJ . „ ,, ,. . ,, . , 

J J ° ° velocity dispersion a (lor further d iscussion m this topic see 

such emission lines we can, to some extent, infer the physi- | M elnick et aUM iTelles and Terlevichlll993h . One funda- 

cal properties inside the regions, the virial mass of the sys- , . , , . , , , ,. ,. . r 4 . . . . 

mental relation reported by these authors is L ex a , which is 
tern, temperatures, densities, and their ages. One important . , ., „ , -, , . .. „ ... , . , , 

° equivalent to the labcr-Jackson relation tor effipticaf galax- 

parameter for deriving dynamical properties is the veloc- . ttt~, — — r^r — i IHT^K n- n .. , • 

, _ TT7TT -_, , — : ics (labor and Jackson lym). Given the potential mipor- 

lty dispersion of the emission line a = FWHM/2v21n2, r n ,. ,. . .. ' , ,. , , T , 

fl ' tance of hnding a valid direct relation between L and a to 

where FWHM is the full width at half maximum of the line. , , ,, . , , , . r TT . , , , . 

understand the mternaf dynamics of H 11 regions but that 
By studying the integrated Ha emission-line profiles of gi- , , , , . , , . , .. ., 

J J ° ° 1 -r- c 1 1 "1 , observations have not seemed to yield a unique relation, it 



ant extra-galactic H 11 regions, ISmith and Weedmanl (|197CJ . 
found velocity dispersions, a ~ 25 km s" 1 . Since then 



has been the subject of some controversy. 



After an initial claim by ITerlevich and Melnickl (|l98ll , 
that they had observed a relation L oc a 4 implying that H 11 
* E-mail: javier@astro.su.se regions are in virial equilibrium, 
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iRov et all (|l986f l and iHippeleinl (|l986l ) found a similar 
relation, but with exponents ~ 3 and ~ 6, respectively, while 
iGallagher and Hunterl (|l983T ) found n o clear proporti o nality 
between L and a. On the other hand. lMelnick et al.l (|l987r i 
adduced problems in the selection of the sample made by 
Gallagher & Hunter, si nce they used regions with lo w ve- 
locity dispersion, while iTerlevich and Melnickl (|l981T ) used 
giant H H regions, with high velocity dispersion [a > 13 
km s~ ) . iMelnick et al.l (|l987l ) reported a new value of ~ 5 
for the exponent in the L — a relation. A common character- 
istic of all these studies is that they use at most few tens of 
the brightest regions, most of them in the high velocity dis- 
persion range , taken from different galaxies, i.e., introducing 
extra uncertainties from the assumed distances to different 
systems. This produces additional scatter as extinction and 
distances are different for each galaxy. 

The obvious next step was then to study all the re- 
gions available from the same galaxy where the distance 
would have no effect on the scatter a nd the extragalac t ic and 
Galactic extinction would be similar. I Arsenault et al.l (jl990|) 
studied 127 regions in the galaxy NGC 4321 and found no 
clear L — a relation, except for the subsample of the highest 
surface brightness regions, for which they report an expo- 
nent of ~ 2.6. Probably the mo st remarkable aspect of the 
study bv lArsenault et al.l (| 19901 ) is that they define an up- 
per envelope for the log(L) — log (a) diagram, an upper limit 
for the luminosity at a certain ve l ocity widt h. This enve- 
lope w as also found bv lRozas et al.l |l998|) and lRelafio et al.l 
(2005) with slopes of 2.6 and 2.0 respectively, and again by 
Rozas ct al. (2006a) who found 2.9 for a sample of super- 
giant H n regions. All these papers relate the upper envelope 
to where the regions in virial equilibrium are located, while 
broader line profiles have contribution from non-equilibrium 
processes 

The line profiles of spatially unresolved H II re- 
gions have been traditi o nally fitted by Gaussians 
JSmith and Weedmanl Il970l ; ITerlevich and Melnickl Il98ll; 
lArsenault et al. | Il990l. and man y others). iRozas et al.l 



(2006gj and IWestmoquette et all (|2007l ) distinguish a 
broad low intensity component added to the main spectral 
component. 

This has been more a question of mathematical sim- 
plicity than of the physical processes which influence the 
line shape. Studies of spa t ially resolved H h regions (e.g 



Meaburn and Walshl Il98ll ; iMeaburnl 1 19841 ; IClavtonl 1 19861 ; 
Chu and Kennicuttlll994l ) showed clearly that not only do 
they receive multiple inputs of energy, but that these in- 
puts vary from place to place (and time to time) within a 
region, so that the resulting integrated emission line profile 
is a luminosity weighted mean of all these inputs. We can 
pick out two principal contributors to the final line shape: 
the turbulent motion within the ionized gas (particularly 
where it interacts with the internal and surrounding neu- 
tral gas), and the emission fr o m dominant expanding shell s 



l|Dvsonlll979l ; IMeaburnl [l980l ; IGallagher and Hunterl 1 1983h . 
The contribution from the turbulent motions is gaussian, 
as it is produced by random motion, while the expanding 
shells, if present, are identified as equidistant in velocity 
from the main peak. These symmetric components produce 
non-gaussian wings or, if the expansion velocity and the lu- 
minosity are high enough, secondary peaks. In addition to 
these two contributions, there may be more local expansion 



due to the winds of multiple in dividual stars and super- 
nova shells JRedman et al.ll2003l ). Also the hyperfine struc- 
ture of Hydrogen split Ha in seven components ~0.14 A 
apart, which should be considered. 

With such a complex ensemble of contributions it is 
importantthe to take into account the shape of the spectral 
instrument response to accurately correct its effect on obser- 
vations. Although, again for the sake of mathematical sim- 
plicity, it has traditionally been assumed Gaussian, the in- 
strumental response is better described by a Lorentzian (e.g. 
lBland-HawthorrJll995l ) . Thus, instead of fitting a Gaussian 
(assuming the t urbulent motion as the dom inant feature in 
the line profile) . iMoiseev and Egorovl (J2008J) concluded that 
Voigt profiles are a more realistic characterization of the in- 
tegrated emission line profiles of H II regions, since it meant 
assuming the instrumental response is Lorentzian, not Gaus- 
sian (more about this will be discussed in sect. 14.2) ). 

During the present study we take particular care of the 
instrumental response and fit Gaussians convolved with it 
to the integrated Ha line profiles, comparing it with the di- 
rect fit to a Gaussian. In sect.[2]we describe the observations, 
sect.[3]explains the data reduction, focusing on the improved 
steps carried out for the present study. The line profiles are 
analyzed in sect. [4] with special emphasis in the L — a re- 
lation. Results and discussion are in sect. [5] and conclusions 
presented in [5] 



2 OBSERVATIONS 

We use a narrow-band Ha image to identify and measure 
the Ha emission flux from the H II regions in M 83, and use 
complementary scanning Fabry-Perot observations over the 
central ~ 2.3 kpc radius to derive the dynamical information 
from their integrated spectra. 



2.1 Ha imaging 

The narrow-band Ha observations were carried out in 
April 2003 with the Danish 1.54 m telescope, equipped with 
DFOSC, at ESO, La Silla, Chile. On the 2kx2k chip, the 
pixel scale was 0'.'39, and the field of view was 13'7xl3.'7. 
We used the narrow-band filters ESO #693 (A e ff « 6562.3 A, 
AA^62.lA) and #697 (A cft ^6654.4 A, AA«61.6A). The 
latter was used for continuum subtraction after matching 
with the Ha image using 17 stars to align the two images. In 
total the exposure times were 2100 s and 1500 s in the two fil- 
ters, respectively. The data reduction was performed in the 
usual manner including bias level correction, flat-fielding, 
background correction, flux calibration, foreground star re- 
moval, and final continuum subtraction. All these reduc- 
tion steps have been performed using various IRAF pack- 
ages. A correction for a Galactic extinction of 0.18 m was 



ap plied. The derived g l obal flux agrees with that obtained 
bv lBell and Kennicuttl (|200ll ) . 

Our Ha filter contains the [N 1 1] lines at 6548 and 
6583 A at either side. The flux ratio of the stronger [N 1 1] line 
(6548 and 6583 A) to that of the Ha line has b een measured 
to 0.45 in the disk of M 83 (|Boissier et alj|2005l ). Considering 
the transmision curve of filter ESO #693 we estimate that 
the [N 1 1] contribution to the measured fluxes in the Ha fil- 
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ter is 25%, and we have corrected the fluxes by a factor of 
0.75 to compensate for this. 



2.2 GHqFaS data 

We obse rved M 83 with the GH qFaS Fabry-Perot interfer- 
ometer IjHernandez et al.l |2008| ) mounted on the Nasmyth 
focus of the 4.2 m William Herschel Telescope (WHT) on 
La Palma. The observations were carried out in two differ- 
ent sets, July 6 of 2007 (~ 60 min) and March 18 of 2009 
(~ 90 min). Using the interference order p=765, GHqFaS 
was tuned to a mean Finesse of ~17 (at Ha) which allow us 
to scan the Ha A6562.78 emission line at 8.2 km s _1 (0.2 A) 
over a free spectral range of 392 km s" 1 (8.6 A). GHoFaS 
was set in low spatial resolution mode (512 x 512 pixels, 
0.4 "/pixel), integrating 5 seconds for each one of the 48 
channels during 38 cycles. The mean seeing was 0.7" ± 0.1" 
during the first run, and undetermined for the second run 
due to a malfunction of the instrumentation of the WHT. 
We do however measure the approximate seeing from our 
final cube image after adding both observations and find a 
mean value of 2 ± 0.2 " across the field. 

The intensity map of the reduced and derotated data 
cube is illustrated in fig.[l] where all the selected H n regions 
are marked (see sect. 13.41 and 14. II for further details). 



3 DATA REDUCTION 

The commonly adopted data reduction for Fabry-Perot 
observations aims at deri ving line-of-sight velocity fields 
over extend ed objects fe.g. lde Vaucouleurs and PenceJll98Cl : 
lAmram et "ail [2002). T he velocity field for M83 was pre- 
sented and analyzed in lFathi et al.l (|2008r ) and here we only 
focus on the shapes of the observed Ha emission line pro- 
files, therefore we need to revise and improve some of the 
data reduction steps to accommodate the accuracy that is 
necessary for our analysis. In the following three subsections, 
we describe the three reduction steps that we have improved: 
derotation, sky subtraction and instrumental response cor- 
rection. 



3.1 Derotation 

As GHaFAS is mounted on the optical table of the Nasmyth 
focus, the Fabry-Perot observations are affected by the rota- 
tion of the field. The optical derotator provided by the Isaac 
Newton Group of Telescopes has the disadvantage of cover- 
ing a field of 2.5' x 2.5', which is almost a factor two smaller 
than the 3.4' x 3.4' field of view of GHoFaS. Using a IPCS 
photomultiplier with rapid readout (|Gach et alj|2002j . |2003| ) 
we are able to store the image from individual channels sep- 
arately. This setup is ideal for derotating the observed cube 
a posteriori without any loss of sensitivity or signal. For our 
observations this means that images with 5 sec interval are 
stored. We rely on the system spee d to record i mages and 
improve the procedure developed bv lFathi et al.l (J2008T) and 
iHernandez et al.l l|2008l ) and apply the following scheme: 

(i) Phase-map correcting all individual channels sepa- 
rately: The phase-map correction rectifies the position of the 
incoming photons for the interference rings in each channel, 



disentangling the wavelength, A, from (x, y) coordinates and 
building one sub-cube (with A as the third coordinate) out 
of each channel. 

(ii) Stacking channel maps: Given the freedom of packing 
together any arbitrary number of consecutive channels, we 
ensure that two criteria are fulfilled: (1) we are able to detect 
at least three bright point sources in every group of channels, 
and (2) we see no sign of field rotation within them (ratio 
between x and y axis for point sources show less than 15% 
difference in the stacked images). An image is built for every 
such group and is a nalyzed w ith FIND procedure (part of 
IDL astrolib library, lLandsma n 1993) until our criteria are 
fulfilled. 

(iii) Searching and storing: The coordinates of point 
sources detected in all the stacked images are calculated by 
determining the center of the point sources using a centroid. 

(iv) Calculating the rotation and translation: For every 
stacked image we ensure that all identified sources coincide 
with the reference source. We use the POWELL routine of 
IDL, to compute the transformation (Ax, Ay, A9) that min- 
imizes the maximum deviation of any of the point sources 
after correcting. Other criteria, such as minimizing the av- 
erage of the deviations can be easily implemented in the 
code. 

(v) Finally, we apply the calculated rotation and trans- 
lation to all the sub-cubes and adding them in the final 
derotated cube. 

Re- measuring the position of the sources in all the dero- 
tated images, we are able to estimate the errors committed 
in the process. These errors mainly arise from the marginal 
uncertainties in the determination of the center of the point 
sources due to the combination of faint fluxes, smearing of 
the data by seeing, and the finite spatial sampling element. 

It is worth noting that the decision about the integra- 
tion time is critical for the derotation. We are stacking the 
channel maps under the assumption that the field rotation 
can be ignored within the stacked images while at the same 
time achieving enough signal to accurately measuring the 
center of the sources. This equilibrium is even more difficult 
to achieve when some of the point-like sources are emission- 
line objects, as they will appear in a limited number of chan- 
nels. Hence, if the reference points are mainly continuum 
sources, short exposure times are better, while if the dero- 
tation relies mainly in emission-line sources longer exposure 
times are preferred, so that the center of such objects will 
be well measured. 

In the case of M 83 we use 14 H n regions and field 
stars evenly distributed throughout the field and 5 second 
exposure time per channel. We estimate the accuracy of the 
derotation to be ~ 1 pixels, with limit cases of ~ 1.5 from 
comparing the position of the sources in the derotated im- 
ages before adding them in the final cube. Our analysis deals 
with integrated line profiles, co-adding 22-600 pixels, so we 
find the errors induced by the derotation negligible. 



3.2 Sky subtraction 

iDaigle et al.l (J2006J) presented an improved sky subtraction 
scheme by measuring the sky in several places around the 
galaxy, fitting a polynomial and interpolating the sky over 
the galaxy to be subtracted from each channel map. This 
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GHocFaS Cube Image for M83 
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Figure 1. Intensity map of M83. Black circles indicate the selected Hll regions, yellow circles the rejected ones. The center of the image 
corresponds to the coordinates RA:13h37m00.8s DEC:-29d51m54s. 



scheme cannot be used for our data set since M83 (~ 12' 
major diameter) is larger than the GHcvFaS field of view. 
Hence, we measured the sky spectrum by observing one cy- 
cle of a blank field 10' away from the galaxy after every three 
cycles on the object. This strategy ensures that we can mea- 
sure sky variations with time and airmass, and it enables us 
to statistically measure the sky spectrum. As the sky field 
is observed in a position slightly offset from the galaxy, the 
structure (if any) of the sky over the field of view could be 
different from that over the galaxy. We opt for not making 
any assumptions about the sky variation over the galaxy 



and its surroundings, and subtract from the galaxy cube a 
mean sky spectrum. 

We calculate the mean s\ and the standard deviation 
A> of each channel in the sky cube. The mean sky profile 
shows a clear peak at ~ 6577 A, in agreement with the 
po sition of an OH doublet (unresolved in our data) according 
to lOsterbrocket all (|l996l ). 

We then subtract s\ ■ n P i X to correct the profile of an 
integrated H II region with n P i X pixels, which according to 
the central limit theorem, should have a standard deviation 
of A = Aa ■ Jn V i X in the case of a random distribution. 
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Figure 2. Every point in this figure is the standard deviation 
of the number of counts in 2000 groups of n p i x pixels, measured 
in two ways. The groups symbolized by triangles are formed by 
pixels randomly selected across the field of view, and they do 
follow the behavior predicted by the central limit theorem (solid 
curve). The asterisks represent groups where pixels have been 
taken in boxes, so they all belong to the same area of the image. 

Figure[2]demonstrates that this is not the case in our data, as 
it shows the study done for groups of pixels ranging from 20 
to 750 in two different ways. The first one is picking random 
pixels within a channel of the sky cube while the second 
one picks a box with the same number of pixels. Doing this 
procedure 2000 times allows us to make statistics about the 
mean value and the standard deviation of the groups. We 
see clearly in fig. [2] that the randomly picked pixels follow 
the expected behavior marked by the central limit theorem, 
while the boxes have bigger scatter due to different areas 
having systematically different number of counts, i.e. due 
to the sky having structure. As we are subtracting the sky 
to single integrated H II regions, the latter scatter provides 
us with a reliable measurement of the error committed in 
subtracting the sky. 

3.3 GHqFaS instrumental response 

We measure the Instrumental Response Curve (IRC) of 
GH«FaS by observing the calibration Neon lamp line at 
6598.95 A. To this line we fit both a Gaussian and a 
Lorentzian, getting a velocity dispersion for the Gaussian 
a = 8.75 km s" 1 (equivalent to FWHM=20.6 km s" 1 ) and 
FWHM=18.6 km s _1 for the Lorentzian. As shown in fig. O 
both functions fit well the core of the IRC, giving similar 
results for the FWHM, but in the wings the Gaussian per- 
forms clearly worse than the L orentzian, as has been known 
for many years (e.g. iBland-Hawthornll 19951 ). 

We will demonstrate in sect. 14.21 that these residual 
wings are indeed very important for calculating the shape 
of the profiles. 



3.4 Building a catalogue of H n regions 

To detect and catalogue the H n regions in M 83, we use 
the software packa ge REGION written by Clayton Heller 
(jRozas et alj|2000l ). This package is applied to the reduced, 
flux calibrated and continuum subtracted Hcv image and 
identifies the regions, measures the position of the center 
of each region, obtains the area of the region in pixels, and 
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Figure 3. Instrumental Response Curve of GHoFaS : asterisks 
represent the calibration line NeA6568.95 A , the dashed line 
is the fit to a Gaussian and the continuous line the Lorentzian 
fit. Although both lines are correct for the core of the line, the 
Lorentzian performs better in the wings. 



calculates the total Ha flux, integrating over all the pixels 
in the region and subtracting the sky value. The part of the 
program that calculates the H II region flux is capable of 
dealing with an irregular background, so we define 43 cir- 
cular areas where the program estimates the average value. 
The background value applied to each H n region is that 
corresponding to the mean value of the nearest sky area. 

The analysis is done in an automatic process, but some 
fine adjustments to the regions are allowed in an interac- 
tive way, intended to modify their size by manual editing. 
So, with the help of the user, the "apertures" can have an 
arbitrary shape and not precisely circular, as seen in fig. [4] 

In order to calculate the boundaries of a region, we cal- 
culate the standard deviation in the nearest background 
and extend the aperture until the surface brightness de- 
creases to a value of twice the standard deviation over 
the background level. Any pixel with a value lower than 
2.5 x W~ 16 ergs~ 1 crn~ 2 is considered part of the diffuse ion- 
ized gas, hence the flux integration within the aperture is 
done only with pixels whose signal is higher than this critical 
level. In addition, all the regions with areas smaller than 20 
pixels or integrated flux smaller than 3.5x 10 -15 ergs' 1 cm~ 2 
are discarded. 

In some cases it is difficult to decide whether a feature 
is the product of a single H II region or the overlap of two 
or more smaller regions. In those cases, a plot of the image 
surface is produced. If more than one peak is visible, with a 
valley deeper than 75% of the height of the smallest peak, 
they are separated and treated independently. Otherwise, 
the reagion is treated as a single region. 

Using these criteria 1029 regions are found over the face 
of the galaxy, of which ~ 350 are inside the region covered by 
the GHoFaS data cube, i.e. the central 2.3 kpc radius. We 
cross match the H n region catalogue with the cube image 
from GHoTaS, and find that 193 regions can be identified 
in the Fabry-Perot cube. We further note that the GHaFAS 
observations are not as deep as the Ha image, and some re- 
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lobs is formed by the intrinsic emission line profile from the 
H n region J rea i convolved with the IRC: 



Figure 4. Example of H II regions selected by the program RE- 
GION. The vast majority of the borders were selected automat- 
ically by the program with the parameters determined after sev- 
eral iterations. In a few cases the borders were edited manually 
especially to separate some overlapped regions. 



gions have to be omitted from our analysis since they cannot 
be identified clearly. 



4 ANALYSIS 

4.1 Extraction of line profiles 

Due to the fact that the Fabry-Perot data acquisition is in- 
trinsically different from the narrow-band image, we have 
restricted our study to those regions which displayed an in- 
tegrated line profile with amplitude-over-noise A/N > 10. 
The A/N ratio has been measured fitting a single Gaussian 
to the line profile, been A/N the amplitude of the Gaussian 
divided by the standard deviation of the spectrum outside 
3<r from the peak. Although we will discuss in the next sec- 
tion that fitting to a Gaussian is not valid when studying the 
shape of an integrated line, it is a satisfactory approximation 
to make a A/N estimation. 

Applying these two selection criteria leaves us with a fi- 
nal sample of 157 H II regions for which we have good quality 
Ha emission line measurements. Eighteen of those regions 
are presented as an example in fig. [5] 

As illustrated in fig. [TJ the selected regions are dis- 
tributed over the entire GHoEaS field, and thus allow us 
to study the very center of the galaxy, the arms, inter-arms, 
isolated H II regions as well as those in groups. 



Jobs — -*n 



IRC, 



(1) 



where I rea i contains the thermal and non-thermal broaden- 
ing. 

To study the non-thermal broadening effect, a cor- 
rection for the contribution of the instrumental and ther- 
mal broadening is needed. It is common to assume 
that the thermal broadening and the IRC are described 
by Gaussians wit h velocity disper si ons cr t h and o-j nst , 
respectively (e.g. iRozas et al] 1 19981 : iRelano et alj 120051 : 
Martfncz-Dclg ado et al] 12007V In such idealised case, the 
convolution of Gaussians leads to quadratic co-addition of 
the velocity dispersions and thus the non-thermal velocity 
dispersion can be calculated as: 



2 2 2 2 2 

Cnt — °obs — °mst — °"th — 0"n, 



(2) 



where a- lns t = 8.75 km s (see sect. I3.3f) and, assum- 
ing an electron temperature of 70 00 K for the H II re- 
gions of M83 (JBresolin et al] I2005J ). a th = 7.63 km s" 1 
(for details about the calculati on of thermal broadening 
see lO'Dell and Townslevl 119881 ). The last term, a n cor- 
rects f or the fine st r uctur e of the Ha line, described in de- 
ta il in Clegg et al] II 19991) and , for very narrow Ha lines, 
in ICumming and Meikld II 19931 ). Simulations based on the 
seven fine-structure components for N e ~ 10 2 and T e ~ 7000 
K yields a value for this correction of a n — 3.2 km s _1 . 

The main caveat of the quadratic subtraction is the as- 
sumption that the IRC of a Fabry-Perot interferometer is 
well represented by a Gaussian, when it is better described 
by a Lorentzian. 

Figure [3] shows that although a Lorentzian is a better 
representation than a Gaussian, there are asymmetries in the 
GHaFAS IRC. Therefor, we prefer not to assume any analyt- 
ical function, using instead the measured non parametrised 
IRC, averaged over the full field of view. We note further- 
emore that the IRC does not vary considerably across the 
field. Only for comparison reasons, all 157 H II regions were 
also fitted according to eq. @, although we stress once again 
that this is a poor approximation. 

The fits are carried out using the IDL procedure MPFIT 
(JMarkwardt 2009), a robust routine for \ 2 minimization us- 
ing the Levenberg-Marquardt technique. When the best fit is 
achieved, the result is four parameters (seven if two Gaus- 
sians are fitted): continuum, central velocity (s), height (s) 
and velocity dispersion(s). Each of these parameters has a 
formal standard deviation around the optimal value, which 
allows us to perform a Monte-Carlo simulation in order to 
estimate the scatter expected around the model in each of 
the channels. 



4.2 Fitting procedure 

We analyse the physical properties of the H II regions by 
fitting single and multiple Gaussian components to the inte- 
grated Ha emission line profiles. An observed emission line 



5 RESULTS AND DISCUSSION 

An example of the fit provided by the procedure described 
in sect. l4.2l is presented in fig. |S]for one and two Gaussians, 
both using direct Gaussian fitting and the convolved model. 
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Figure 5. Examples of Hll region integrated emission line profiles extracted as described in sect. 14.11 Top three files are the 9 regions 
with worst A/N, while bottom three rows show the 9 regions with best A/N. 
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Figure 6. Example of the result of the fitting procedure (sect. [42}- Upper left: Fit to one Gaussian (blue) and one Gaussian convolved 
with the instrumental response (red). Center left: Fit to two Gaussians, the black line being the addition of both components. Lower left: 
Fit to two Gaussians convolved with the instrumental response. Again the black line is the addition of both components. Right panels: 
Residuals to the fits. Vertical black lines are the expected scatter in each channel due to the addition of three sources: scatter in fitted 
parameters, sky subtraction and Poisson noise of the incoming photons. 



5.1 Convolution versus quadratic subtraction 

The reduced y 2 of the fits are ranging from 0.8 to 10 for both 
methods, in 76% of the regions favorable to the convolution 
method. The direct Gaussian failed to fit the two most lu- 
minous points in 65% of the profiles due to the weight that 
the wings have in the fitting, while the convolved Gaussian 
give a better fit in both the wings and the core. 

On the other hand, reduced \ 2 is U P to 10 because our 
noise is underestimated, as shown by the fact that the dis- 
persion in the continuum is bigger than the individual un- 
certainties of any of the points. When we use the dispersion 
of the luminosity of the points in the continuum (> 3cr away 
from the center of the profile) as the noise for these points, a 



X 

while the convolved model reduces to the range 0.7-2. Thus, 
we can conclude that for the case of one Gaussian, the fit 
using the convolved method is more robust. 

When performing a two-Gaussian fit to our data, the \ 2 
of the fit is improved, but also the uncertainties in the de- 
rived parameters are bigger. The formal uncertainties given 
by MPFIT when fitting a single Gaussian are 1-3% for the 
height and 2-8% for the velocity dispersion in the case of 
single Gaussian fitting, while for the case of the convolved 
Gaussian the uncertainties are 4-8% and 2-12%, respec- 
tively. When fitting two Gaussians the uncertainties increase 
to ~ 40%, even larger when the convolved model is used. 
This is explained by the fact that the line profile is present, 
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Figure 7. The velocity dispersions arc displayed against the lumi- 
nosity for both methods. In red the velocity dispersion calculated 
by direct fitting of one Gaussian, once corrected from thermal and 
instrumental broadening. In blue, the velocity dispersion of the 
same regions calculated with the convolution of a Gaussian with 
the instrumental response. On the bottom panel the subtraction 
of both sigmas, giving values around 8 km s for all the regions. 



in the best of the cases, in a dozen channels and we are 
fitting six parameters to it. 

When fitting two Gaussians requiring amplitude ratio 
greater than 1 : 5, 32 regions can be fitted with less than 
20% error in the velocity dispersion, and when two convolved 
Gaussians are fitted, 11 regions can be fitted following the 
same criteria. Thus, we conclude that our observed line pro- 
files are adequately fitted by one single Gaussian only, and 
inferring an additional Gaussian, will result in very weak 
parameter determination for both components. 



This relation has been a matter of discussion for the 
last four decades. Table [l] summarizes the main articles re- 
lated with the L — a relation, the sample sizes and the re- 
sults they reported. All these studies are based on Fabry- 
Perot observations, and used the quadratic subtraction to 
correct from the instrumental response. While the ear- 
lier works pointed to a correlation between both variables, 
with slope between 2 and 7, as observations were refined 
studies showed that the best linear relation found when 
plotting log(cr) ag ainst log(-L) is that to the upper enve- 
lope of the p oints (Arsc nault et al.lll990l ; lRozas et al.lll998l : 
IRelano et a l. 2005). These latter works showed a scatter di- 
agram very similar to fig. [8j where we present the L — a rela- 
tion for both the convolved and the quadratically subtracted 
method, assuming for the Luminos i ty a d istance of 4.5 Mpc 
(|Thim et alJl2003h . lBeckman et al.l (|2000h described this en- 
velope as the locus of the regions with minimum velocity 
dispersion for a given luminosity. Such a minimum velocity 
dispersion suggested that those H n regions near the enve- 
lope are in virial equilibrium, which implies that the great 
majority of the regions are not in virial equilibrium. This 
is entirely plausible considering the many contributions of 
dynamically impelled gas to the emission line profiles in gen- 
eral, as we have seen in s e ct. [T] 

As lArsenault et all (J1990T I. iRozas et all Jl99ct ) and 
IRelano et al.l l|2005l ). we do find an upper envelope for the 
log(L) — log(cr) diagram. Our data do not allow us to dis- 
cuss whether virial equilibrium can be applied to the re- 
gions near the envelope, but only to report that such an 
envelope, with minimum velocity dispersion for a given Lu- 
minosity exists, as claimed by the works cited earlier. The 
slope of this envelope depends on the method used to cor- 
rect the influence of the IRC. For the convolved model we 
find log(L Ha ) = (0.9 ± 0.2)log(a) + (38.1 ± 0.1), and we 
insist on this as the best of the two methods, as it takes 
into account the full shape of the IRC. On the other hand, 
and to be able to compare with previous studies, the enve- 
lope for the quadratic subtracted a is also calculated, giving 
log(L Ha ) = (2.0±0.3)log(a) + (36.4 ±0.4), in perfect agree- 
ment with the latest results in Table [T] 



5.2 Fitted parameters 

Comparison of the fitted parameters for both methods show 
that the position of the center of the line is invariant, an 
expected result that confirms that in order to get velocity 
maps it is enough to assume a Gaussian for the line profile. 
Once applied eq.[2]to correct the thermal and instrumental 
broadening from the Gaussian fit, and the thermal broaden- 
ing from the convolved Gaussian, we obtain the velocity dis- 
persions presented in fig. [7] Quadratic subtraction systemat- 
ically overestimates the velocity dispers ion, in full agreement 
with the recent study presented by iMoiseev and Egorovl 
(2008), which predicted a constant overestimation. In the 
case of GHcvFaS and with the setup used here, the differ- 
ence is ~ 8 km s" 1 (cf., fig. [Tf. 



5.3 L — a relation 

The relation between luminosity and velocity dispersion, can 
be described by: 



log(L) = a ■ log(o-) + b 
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(3) 



5.4 Hn regions as a possible distance estimator 

Since the earliest observational studies of the L — a relation 
the idea of using it to help formulate a d istance indicator has 
been discussed by severa l authors (e.g. iMelnick et aLlll987l ; 
lArsenault and Rov|ll988f). 



In their article lBeckman et al.l (|2000r i showed that when 
histograms of the number of H n regions are plotted vs. lu- 
minosity, this luminosity function for a given galaxy has a 
break at log(L) = 38.6 (±0.1)rie:r, a break first reported by 
iKennicutt et alj (119891') and repeatedly confirmed by other 



iy i 

authors (e.g. Rand 1992J^ Rozas et al.1 1 199a ; iBradlev et al.1 
120061 ). iBeckman et al.l (|2000T ) put forward the luminosity of 
this break as a possible standard candle, by analogy with, 
but more powerful than the standard candle offered by the 
knee of the luminosity function of planetary nebulae. All 
the above mentioned studies are based on reasonably face- 
on galaxies, ( i < 50°), in order to optimize the number of 
regions detected, prevent crowding effects, and minimize the 
extinction by dust in the disks of the galaxies. These effects 
will limit, to some extent, the utility of the 38.6 dex break, 
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Table 1. Summary of papers relevant for sect. 15.31 All are based on Fabry-Perot data, with quadratic subtraction of the 
IRC (details in sect. l4~2t 



Paper 


Sample 


Results 


Smith and Weedman (1970) 


9 


High velocity dispersions (o > 25km s — 1 ) 


Melnick (1977) 


16 


Rccct 227 


Terlevich and Melnick (1981) a 


12-15 

17 


L(H/3) oc o 4 


Rov et al. (1986) 6 


log(L Ha ) = (3.15 ± 0.83) logo- + (35.37 ± 1.17) 


HiDDelein (1986) 


43 


log(L HQ ) = (6.6 ± 0.4) logo- + (31.1 ± 0.3) 


Melnick et al. (1987) 


22 


log(i H /3) = (4.92 ± 0.30) logo + (33.24 ± 0.31) 


Arsenault et al. (1990 ~[£ 


127 


log(L Ha ) = 2.18 logo + 36.5 


Rozas et al. (1998) C 


200 


log(i HQ ) = 2.61ogo + 36.15 


Relafio et al. (2005) C 


373 


log(L H a) = (2.0 ± 0.5) logo + (36.8 ± 0.6) 



Uses data from Melnick (1977) 

No mention about the IRC correction 

Defines an upper envelope for the L — o relation, not a fit 



notably for galaxies with higher inclinations, but the validity 
of the result for low inclination galaxies appears solid. 

The data displayed in fig. [8] show that the linear upper 
envelope referred to above sets in at log(L) between 38.5 
and 38.6. This value is based on a distance to M 83 based 
on measured Cepheid variables. This is, to our knowledge, 
the first study that shows a break point in the envelope of 
the log(L) — log(cr) diagram. Without venturing into the un- 
derlying physical cause of this change in slope its existence 
suggests that there is such a cause, and tends to reinforce 
the use of the log(L) — log(N) histogram, or indeed the enve- 
lope of the log(L) — log(o-) itself, with their respective break 
points, as possible versions of a secondary standard candle. 

The uncertainty in the determination of the point where 
the change of regime happens in our data yields errors of 
the order of 0.4 Mpc, compar able with the 0.2 M pc error of 
the Cepheid measurement by iThim et a.1.1 |2003). A better 
determination of this point for M 83 should be possible, given 
the fact that our 3.4' x 3.4' GHaFAS field of view covers only 
~8% of the disk in M 83. 



5.5 Ha equivalent widths 

Following the formalism developed in the Starburst99 
(SB99) model (jLeitherer et a.1.1 Il999f ) some crude assump- 
tions can be used to obtain a rough idea about the ages of 
the observed HII regions from our data. We acknowledge 
that the situation is very complex not only due to obser- 
vational difficulties but also due to dust scattering effects 
on the continuum levels, nevertheless, as a mere exercise 
we can assume comparable dust content in the spiral arms 
and the bar, an instantaneous burst of star formation for 
the HII regions, and a constant a verage metallicity, Z = 
0.04 IjBresolin and KennicuttJ 120021 ) and use the measured 
Ha equivalent widths to estimate integrated average ages 
of H n regions in M 83. Comparing the SB99 models with a 
typical value of the EW of 60 A for those HII regions re- 
siding in the spiral arms, we obtain a mean age of 9 Myr. 
A similar exercise for the H II regions in the bar yields an 
average EW of 20A , which corresponds to a mean age of 12 
Myr. Although these values are lower limits, the systematic 
difference between the EWs in the arm and bar region sug- 



gests that the HII regions in the spiral arms are typically 
younger than those in the bar of M 83. Leakage from other 
orders of interference, due to imperfectly efficient interfer- 
ence, affect both our measured Ha continuum luminosity, 
resulting in an underestimation of the EWs since the equiv- 
alent width (the ratio of the Ha luminosity to that of the 
nearby continuum) is very sensitive to the continuum value, 
Fabry-Perot observations should produce values for the EW 
which are systematically low, so the actual ages should be 
interpreted with caution. 



6 CONCLUSIONS 

We have presented an improved method for analysing emis- 
sion line profiles from Fabry-Perot interferometers, and have 
applied it to GHaFAS Fabry-Perot data for the central 2.3 
kpc radius of the nearby Sc galaxy M 83. We have considered 
the true instrumental profile without any parametrisation, 
and derived observed emission line profile parameters by a 
convolution of a Gaussian with the instrumental response. 

From a narrow-band Ha image, a catalog of H II regions 
was extracted and we were able to determine the flux, posi- 
tion, and area for 1029 star forming H II regions. Improving 
some of the Fabry-Perot data reduction procedures (sect.[3| 
we have measured accurate integrated Ha emission line pro- 
files for 157 H ii regions in the GHaFAS field. 

Convolving with the IRC before fitting is the more re- 
alistic approach and clearly gives better fitting results in 
terms of smaller \ 2 f° r 76% of the regions. When fitting 
two Gaussians, the errors in the determination of the line 
parameters increase to more than 40%, due to low num- 
ber of degrees of freedom, and for this reason we ruled 
out successful multiple-Gaussian fitting. Following the com- 
monly used procedure, when fitting Gaussians in a direct 
way, the IRC was assumed a Gaussian and quadratically 
subtracted. The direct fitting method leads to an overesti- 
mation of a by a constant value of ~ 8 km s" 1 . We used 
the results from our improved line fitting approach to study 
the L — a relation and found no significant linear correlation 
between these two variables, but an upper envelope defined 
by log(Lna) — 1-3 • log(er) + 37.6. Our data is fully consis- 
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Figure 9. Ha luminosity and Equivalent width map of all the 157 regions within the GHoFaS field of view. 



tent with the previ ous e stimations of Arsenault et al.l (|l990T) , 
iRozas et al.l (|1998t ) and lRelafio et al.l (2005|), when applying 
direct Gaussian fitting. 

The upper envelope emerges in the log(L) — log(cr) di- 
agram at luminosities higher than W 385 erg s^ 1 . 

iBeckman et al.l (|2000h suggested that this luminos- 
ity might correspond to the transition from an ionization 
bounded regime to a density bounded regime (i.e. from 
optically thick to optically thin). We now understand, in 
part from work bv lGiammanco et al.l i|2004r ) that the simple 
Stromgren sphere scenario on which this concept is based is 
too simple, and that the strongly inhomogeneous structure 
of H II regions must be taken into account when explaining 



any product of radiative processes within them. Nevertheless 
an explanation in which there is an increase in the escape 
fraction of ionizing photons from H II regions above this lu- 
minosity is still a reasonable way to account for a variety 
of observations, even though more sophisticated modelling 
will be needed to give final quantitative answers. Our sam- 
ple has extended previous H II region analyses, by including 
a significant number of faint H II regions with their corre- 
sponding values of a, and this has allowed us to better pin 
down the luminosity at which a clean upper envelope in the 
L — a relation emerges. Using the measured value as a sec- 
ondary standard candle we find that the uncertainty in the 
position of the break point yields a contribution to the un- 
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Figure 8. The L — a relation for both methods: convolved and 
direct Gaussian fit. Each point represents one region, the solid 
line is the fit of the upper luminosity envelope. A clearly different 
trend is found for the regions with lower velocity dispersion, as 
showed by the dashed line. Note the difference in the x-axis for 
both panels. 



certainty in the distance to M 83 (0.4 Mpc) which is not 
much greater than the uncertainty in the Cepheid distance 
itself (0.2 Mpc). As our observations cover less than 10% of 
the optical disk of M 83 this suggests that a more complete 
data set could be used to reduce this uncertainty. 
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